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A recently introduced particle-based model for fluid flow, called Stochastic Rotation Dynamics, 
can be made Galilean invariant by introducing a random shift of the computational grid before 
£1h collisions. In this paper, it is shown how the Green-Kubo relations derived previously can be 

resummed to obtain exact expressions for the collisional contributions to the transport coefficients. 
It is also shown that the collisional contribution to the microscopic stress tensor is not symmetric, and 
. that this leads to an additional viscosity. The resulting identification of the transport coefficients for 

^vj the hydrodynamic modes is discussed in detail, and it is shown that this does not impose restrictions 

on the applicability of the model. The collisional contribution to the thermal conductivity, which 
becomes important for small mean free path and small average particle number per cell, is also 
derived. 
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I. INTRODUCTION 

Particle-based simulation techniques have recently become an attractive alternative to more traditional methods for 
the coarse-grained modeling of a fluctuating solvent. A particularly appealing algorithm, introduced by Malevanets 
and Kapral EL 13, and later called multi-particle collision dynamics 0, 0, E| or Stochastic Rotation Dynamics 
(SRD) 0, Isl 19. llOL Hit Il2l IT3I ll4T | provides a "hydrodynamic heat bath," the details of which are not resolved, 
but which provides the correct hydrodynamic interactions and thermal fluctuations. The coarse-grained dynamics 
■^J- ■ ignores molecular details, but enables simulations which span much longer time scales than can be addressed u sing 
molecular dynamics techniques. It has been used to stud y th e behavior of polyme rs JM [Til ITsj . colloids 0, @, Il9| | 
(including sedimentation [lj, [l5j), vesicles in shear flow |2(|, and complex fluids |2lJ [23. In addition to SRD's 
numerical advantages, its simplicity has made it possible to obtain analytic expressions for the transport coefficients 
which are valid for both large and small mean free paths, something which is very difficult to do for other mesoscale 
particle-based algorithms. 

In its original form 0,0 j the SRD algorithm was not Galilean invariant at low temperatures, where the mean free 
path, A, is smaller than the cell size a. However, Galilean invariance can be restored by introducing a random shift 
0, H| of the computational grid before every multi-particle interaction. In addition to restoring Galilean invariance, 
1 ' this grid shifting procedure accelerates momentum transfer between cells and leads to a collisional contribution to 
the transport coefficients. Two approaches have been used to analyze the resulting algorithm and determine the 
shear viscosity and thermal diffusivity. In Refs. ^2] an d a non-equilibrium kinetic approach is used to derive the 
transport coefficients. In Refs. a discrete-time projection operator technique was utilized to obtain Green-Kubo 

relations |23l I24I |25| for the model's transport coefficients, and explicit expressions for the transport coefficients were 
derived in accompanying papers pi, [ToL Hi) . The two approaches are complementary and, for the most part, agree 
in their conclusions. The first is rather straightforward and intuitively appealing, but makes several assumptions 
which are not easily verified. The current approach justifies in detail several assumptions used in the non-equilibrium 
calculations of Refs. ^3 and 0]; it can also be used to analyze the transport coefficients of the longitudinal modes, 
namely the bulk viscosity and thermal diffusivity, which are hard to calculate analytically in the non-equilibrium 
approach Note, in particular, that the collisional contribution to the thermal conductivity has not yet been 

determined using the non-equilibrium methods. 

In this paper, we show in detail how the time series in the Green-Kubo relations for the transport coefficients can 
be resummed in such a way as to eliminate all dependence on the space- fixed cell coordinates of the particles pfj). 
This leads to a dramatic simplification of the sums, and enables an exact evaluation of the collisional contribution 
to the transport coefficients. It is also shown that there are only pure kinetic and collision contributions to the 
transport coefficients, as was implicitly assumed in the calculations of Refs. 0,0]. Explicit expressions for the 
collisional contributions to the viscosities and the thermal diffusivity are given, and the consequences of the fact that 
the collisional contribution to the microscopic stress tensor is not symmetric are discussed in detail. In particular, 
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it is shown that this lack of symmetry leads to a slight change in the longitudinal viscous transport coefficient. The 
hydrodynamic modes are, however, not affected, and it does not impact on the applicability and stability of the 
method, even at low temperature. 

II. SRD MODEL 

In the SRD algorithm, the fluid is modeled by particles with continuous spatial coordinates rj(i) and velocities 
Vj(t). The system is coarse-grained into the cells of a regular lattice with no restriction on the number of particles in a 
cell. The evolution of the system consists of two steps: streaming and collision. In the streaming step, the coordinate 
of each particle is incremented by its displacement during the time step, r. Collisions are modeled by a simultaneous 
stochastic rotation of the relative velocities of every particle in each cell. As discussed in Refs. and §, a random 
shift of the particle coordinates before the collision step is required to ensure Galilean invariance. All particles are 
shifted by the same random vector with components in the interval [—a/2, a/2] before the collision step. Particles are 
then shifted back to their original positions after the collision. If we denote the cell coordinate of the shifted particle 
* by the algorithm is summarized in the equations 

r<(t + r) = r,(t)+Tv,(i) (1) 
v,(f + r) = u[£(* + t)} + + t)} ■ {v t (t) - ufcj(t + r)]}, (2) 

where denotes a stochastic rotation matrix, and u(£?) = jj X^es s Vfc ^ s ^ ne mean velocity of the particles in 

cell £ s . All particles in the cell are subject to the same rotation, but the rotations in different cells are statistically 
independent. There is a great deal of freedom in how the rotation step is implemented, and any stochastic rotation 
matrix consistent with detailed balance can be used. In two dimensions, the stochastic rotation matrix, uj, is typically 
taken to be a rotation by an angle ±a, with probability 1/2 (see Refs. 0,H,@)- I n three dimensions, two collision 
rules have been considered. In the first, (Model A in Ref. [lflj). one performs rotations by an angle a about a randomly 
chosen direction, where all orientations of the random axis occur with equal probability. In the second, (Model B in 
Ref. |lC|). rotations are performed about one of three orthogonal rotation axes, i.e. x-, y- and z-axes of a cartesian 
coordinate system. At each collision step one of these three axes is chosen at random, and a rotation by an angle ±a 
is then performed, where the sign is chosen at random. 

III. TRANSPORT COEFFICIENTS 

Because of the cell structure introduced to define coarse-grained collisions, angular momentum is not conserved in 
a collision 28, 29]. As a consequence, the macroscopic viscous stress tensor is not, in general, a symmetric function 
of the derivatives V a p = <9 Q v^. Its general form can be determined as follows. Both the macroscopic viscous stress 
tensor, a a p, and the velocity gradient tensor, V a p which appear in the Navier-Stokes equation are rank two tensors. If 
the velocity gradients are small, we can assume — as is generally done — that the momentum transfer due to viscosity 
depends only on the first derivatives of the velocity, so that 

cr a p = CaPipVyp. (3) 

Symmetry arguments can be used to reduce the number of independent elements in the compliance tensor C a /3j P - 
Because of the simple cubic grid structure used in the algorithm, we have at least cubic symmetry. In this case, it 
can be shown 29] that the most general form for C a p~ ip is 

C a p 1P l ' p = a 5 a p5 lp + b 5 ai 5(3 P + c SapSp^ + eY a p lpi (4) 

where Y is the rank four unit tensor p9|. a, b, c and e are viscosity coefficients, and p is the mass density of the fluid. 
It follows that 

Vap/p = a 8 a pd\v\ + b d a vp + c dpv a + er Q(37P <9 7 t> p 

= vi(d a vp + dpVa - ^<S Q| 3<9a«a) + v 2 (d/3V a - d a vp) +^5 a pd\V\ + eT a p yp d~ / v p , (5) 

with kinematic shear viscosities V\ = (b + c)/2, v 2 = (c — b)/2, and bulk viscosity 7 — a + (b + c)/d, where d is the 
spatial dimension. v 2 is the viscous transport coefficient associated with the non-symmetric part of the stress tensor, 
and e is a viscosity coefficient related to a possible lack of full rotational symmetry; both v 2 and e are usually zero 
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in simple fluids. If we define a new tensor from the diagonal elements of V a p, namely R a p = T a p lp V lp , the resulting 
form of the momentum equation for a fluid reads 



9v 

^ + (vV)v 



-Vp + V ■ (7 



-Vjj + pOi + u 2 ) Av + p ( (1 - - z/ 2 + 7 ) V(V ■ v) + eV ■ R 



(6) 
(7) 



where p is the pressure. It can be seen from Eq. J7J that effective shear viscosity is v = v\ + v 2 . In momentum space, 
the resulting linearized Navier-Stokes equation can be written as 



d t v a (k) 



d a p 



Aa/3(^i,^2,7,e; k )^( k ): 



where 



A ttj g(i/i, u 2) 7, e; k) = v x {& a fi + - d ^ k a k fJ ^j + v 2 (s a/3 - k a k^ + jk a k + e k^k p T a p lp . 



(8) 



(9) 



A a /3(^i, ^2, 7, e; k) is the matrix of viscous transport coefficients. In a simple liquid, e = (because of invariance with 
respect to infinitesimal rotations), v — i/u and v 2 = (because the stress tensor is symmetric in d a vp). In this case, 
Eq. @ reduces to the well known form [8| 



A,, ;{i> i . ;/ 2 . ~ . < : k! = /' ( 8 a [3 + < ^—^ L k a kf } ) + -/,-,,/•' ,. 



As shown in Ref. @, the discrete Green-Kubo (GK) relation 

oo ' 

A Q/ 3(^i,^ 2 ,7 ; e; k ) = jYfc^y (^ o 'aA(0)|fcA'Cr/3A'( nT )) 



(10) 



(11) 



for the SRD model can be used to expresses the matrix of viscous transport coefficients in terms of a sum of time 
correlation functions of the reduced fluxes I\+ a (k, t) 



k x o aX {t) = 7i +Q (k, t) = -J2( ~ha(t)k ■ A^(t) + Av ja (t)k • A£*(i)] 



rk, 



(12) 



for a = 1, . . . , d, with A^- (nr) = ^ ([n + 1] r) - ^ (nr), A^ (nr) = ^ ([n + 1] r) - g ([n + 1] r), and Au^nr) = 
v x j {[n + 1] t) — v x j(nr). ^(nr) is the cell coordinate of particle j at time ?ir, while £j is it's cell coordinate in the 
(stochastically) shifted frame. The prime on the sum indicates that the t = term has the relative weight 1/2. The 
sum in Eq. I|12|) runs over all N particles of the system. 

The corresponding expression for the thermal diffusivity is Q 



D T = 



c p Nk B T 2 



^2 (kxX\(0)\k\'X\'{n.T)), 



with the energy flux 

XA(nr) = ~J2 



(c v T - u](nr)/2) A£ iA (nr) - -A^ 2 ( nr )A£j A (nT) + rfc B T^ A (nr) 



(13) 



(14) 



where Ct, = dfcs/2 is the specific heat per particle at constant volume of an ideal gas and Aw|(nr) = Vj ([n + 1] r) — 
u|(nr). The thermal conductivity, k, is related to £>t by k = pc p Dx- Here and in the following we have set the 
particle mass equal to one. 
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IV. RESUMMED GREEN-KUBO RELATIONS 



The straightforward evaluation of the GK relations 23, 24, 25] presented in Ref. [9J leads to three contributions 
to the transport coefficients, which were called the kinetic, rotational (or collisional), and mixed terms. The term 
"rotational" and the superscript "rot" refer to contributions from the collisions — stochastic rotations of the relative 
particle velocities — in the SRD-model. For large mean free path, A ^> a, A = Ty/ksT, the assumption of molecular 
chaos is valid, and the kinetic contribution could be determined explicitly. For mean free paths A smaller than the 
cell size a, however, there were finite cell size corrections, and it was not possible to sum these contributions in a 
controlled fashion. The origin of the problem was the explicit appearance of A£ in the stress correlation functions. 

In fact, the appearance of A£ is troubling, since one would not expect this if the cell shifting procedure really does 
restore Galilean invariance. The key to resolving this dilemma is to realize that a proper resummation of the GK 
relations removes this dependence. Consider first the time series X)^°=o k x a aX (nT). By canceling ^-dependent terms 
in successive contributions to this series, it can be shown that 

-k x [A aX (0) + A aX (r)] + h<T a x(nr), (15) 

n=0 n=a 

where a aX {nr) = a^inr) + a™ x (nr) with 

v*\{nr) = -^[v ja {nT)v jX {riT) - S aX v?{t)/d], ff^f (nr) = -- ^ v ja (n T )B jX (nT), (16) 

3 3 

A aX (T) = ±j2 v ^ A tU°)> ( 17 ) 

3 



and 



B jX (nr) ee ([« + l]r) - £j>r) - Tv jX (n T ) = A^ x (nr) - A^(nr) + A^ s A ([n - l]r) - rv jX (nr). (18) 



Bj X is a new stochastic variable which has very simple temporal correlations describing the geometrical properties of 
the underlying lattice. 

Similarly, it can be shown that 

a aX (0) = A aX (0) - A aX (r) + a aX (0) . (19) 

Using these results, the Green-Kubo relation 1)11(1 for the viscous transport coefficients can be written as the sum of 
two terms, 

^0 ^ \ [A aX (p)Af, x , (0) - A aX {r)A [jXI (r)] + \ [A aX (Q)A px , (r) - A aX (r)A px , (0)] 

1 \ 

+~[A px ,(O)+A0 X ,{T)}a aX (O) + [A aX (O)-A aX (T)}J2 *M nT )) ( 20 ) 



71=0 



and 



Jj-fTY E (*A*«A(0)|*A'*i9A'(nr)) . (21) 

^ n=0 

Stationarity implies that the first term in Eq. I|20(l equals zero, and that the last term reduces to 

oo ' - 

[A aX (Q] - A aX (r)] ^A'(nr) = --[A aX (r) + A aX (O)]a 0X ,(O). (22) 

n=0 

Stationarity and time reversal invariance imply that the remaining term also vanishes. Alternatively, the explicit form 
of A aX , Eq (|17fl . can be used to show that this term vanishes. The expression in Eq. I|2U|I is therefore zero, so that 
the Green-Kubo relation for the viscous transport coefficients is still given by dJ, but with the stress tensor a aX . 



5 



A similar calculation shows that the thermal diffusivity is given by l|13fl . with x*i nT ) replaced by x\( nT ) = 
X k x m (nT)+ X r x ot (nT),with 

X k x m (nr) = ^[-^|(nr)^ A (nr) + k B Tv jX (nr)], xT(nr) = i- £ «*(nr)i? iA (nr). (23) 

3 3 

Note that the new stress tensors do not depend on £, the space-fixed cell coordinates of the particles. 

An alternative way to derive these results is to note that time reversal invariance can be used to rewrite and 
(f 1 31) as sums from — oo to +oo. In this way, the discussion in the preceeding paragraph of the n = term can be 
avoided. 



V. CORRELATIONS INVOLVING B'S 



Bi a {tit) is the a-component of the difference between the change in the shifted cell coordinates during one streaming 
step and the actual distance traveled, TVi a . It has a number of important properties which simplify the calculation of 
the transport coefficients. In particular, it will be shown that all stress-stress correlation functions involving one B in 
the GK relations for the transport coefficients are zero, so that, for example, A a p (v\ ,V2,J, e; k) = A^j™ (v\ ,^2,7, e; k) + 
A^Vi^T^k), with 

OO ' 

A^> 1)2 , 2 , 7 ,e;k) ee £ (k x a k S(0)\h'*$(nT)) (24) 

and 

oo ' 

A^K^ 2 , 7 , e ;k) = {ha r a f{0)\ky^f,(nr)). (25) 

n— 

The kinetic contributions to the viscosity were calculated previously in both 2d and 3d and will not be 
discussed here. Properties of the B-correlations which enable an explicit evaluation of expression l|25f) are derived in 
the following subsections. 



A. Factorization of B — v correlations 



The first one, {Bi a (nr)) = for arbitrary n, implies that on average, the distance traveled by a particle during one 
time step is the average of difference of the shifted cell coordinates before and after the streaming step. This can be 
shown as follows. Consider 

(B<x(0)) = <aW-^(0)-T« te (0)) . (26) 

The ensemble average includes averaging over the initial coordinates and velocities of all particles, as well as averages 
over the shift and collision matrix at each time step. Without loss of generality, assume that at t = 0, the x-coordinate 
of particle i is in the interval [0, a). For na < X = x;(0) +rv ix (0) < (n+ l)a, the average of ^ x (r) = ^(r) — A£?(0) 
over the random shift S at time r, denoted by ( }s T , at fixed particle coordinate and velocity, is (see Eq. (31) of |!j) 

(aWkk = na - (Aa x \x)s T =X- a/2, (27) 

so that 

(B tx (0)) = (-a/2 + x l (0)-a x (0)) ■ (28) 

Finally, averaging over the shift at t = gives (^" x (0)\x)s o = —a/2 + 2^(0), so that (B ix (0)) = 0. 

Similar arguments can be used to show that the cross terms in the GK expressions for the transport co- 
efficients involving one B are zero. For the shear viscosity, these terms involve correlations of the form 
(vj X (nT)vj y (riT)vi X (mT)Biy(mT)}. Consider first the case n = m = 0. Performing the average over the shift 5 
at t = t, the average reduces to 



(^(0)^,(0)^(0)^(0) - a /2 - £?,(0)]). 



(29) 



6 



The average over the shift at t = gives zero because it doesn't affect the particle's initial velocities or positions. 
Consider now m = and n = 1. In this case, perform first the average over the random shift at t — 0. The result is 

(v JX (r)v jy (r)v tx (0m y (T) + a/2 - w (r)]) . (30) 

If the probability of any given configuration at t = t in a shifted cell containing particle i is independent of S, the 
average over the shift at t = r factorizes. This is, in fact, the case since the average in l|30[) entails an integration 
over the initial particle coordinates and velocities at t = 0. In this case, the average over 5 can be performed; since 
(£iy( T ) + a/2 — yi(r))\s vanishes, the result of this averaging is zero. An alternative, more detailed discussion of this 
proof is given in the Appendix. 

The argument for general m and n is similar. Analogous reasoning can be used to show that correlations such as 

(v ix (0)vj X (nT)B iy (0)Bj y (nT)) = (v ix (0)vj X (nT))(Biy(0)Bjy(nT)) (31) 

factorize for arbitrary n. 

B. Autocorrelation of B's 

It is straightforward to evaluate equal-time correlation functions of the B variables. Using the results derived in 
the Appendix and Eq. (36) of ||J, one has 

(Bf x (0)) = a 2 /3 and (B ix (0)B jx (0)) = a 2 /6, for i £ j. (32) 

Correlation functions such as (Bi X (0)Bi X (r)) can be evaluated as follows. Take < Xi(0) < a, A£i X = A£j x (0) = ma, 
and A^ x (t) — na. Here and in the following expressions the argument (0) will be omitted for clarity. Averaging over 
the random shift 5 at t = 2t, 

(B zx (0)B ix (t)) = -([(m+l/2)a-Xi(T)}[A£ ix +AZ s ix (-T)-Aa x -TV ix }) 

+ (A&[A£ to + A&(-T)-A&-™ ia! ]>. (33) 

Since (m + l/2)a — Xj(r) does not depend on random shifts at time and r, the average over Sq and 5 T in the first 
term on the right hand side of Eq. 1|33|) at fixed Xi(0) and Vi x (0) vanishes, so that 

(B ix (0)B ix (r)) - -<(A^) 2 ) + (A&[A£« + A&(-r) - rv ix ]). (34) 

Averaging the second term on the right hand side of this equation over So and <5 r , one finds 

^ pa 00 p[(m-\-l)a — Xi]/T 

(A£f x [A& x + A^(-r) - tv ix \) = - j dx, ^ [(m + l/2)a - x> - Tv lx fw{v lx )dv ix . (35) 

a JO m= _ OQ J(ma- X% )/T 

Comparing with Eqs. (18) and (32) of 9], it can be shown that is equal to -((A£ ix ) 2 )+2(A£ ix A£| i K )+A 2 +a 2 /12 = 
a 2 / 12, where the last equality follows from 

<A& x A^> = i[<A£i)-A 2 ] (36) 

given as Eq. (36) in 9]. 
Finally, using JSBJl, 

(B ix (Q)B ix ( T )) =-a 2 /6. (37) 

The average (Bi x (0)Bj X (r)) can be evaluated in a similar fashion. Take m^a < Xi(0) < (mo + l)a, noa < Xj(0) < 
(no + l)a, A£j X (0) = m^, A£j x (0) = nia, A^ i:E (r) = m 2 a and A.£j x (r) = n 2 a. Averaging over S 2t , one has 

(B ix (0)B jx (r)) = -([(n + n x + l/2)o - ^ (r)] [A£ ix + A&(-t) - A& - ™ <x ]) 

+ (A^[A^ x +A^(-r)-A^-r^]). (38) 

Again, since (no + n% + l/2)a — Xi(r) does not depend on random shifts at time and r, the average over <5o and <5 T 
in the first term on the right hand side of Eq. I|38|l vanishes, so that 

(B ix (0)B jx (r)) = + (A^[A& X + A&(-r) - (39) 
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Using Eqs. (J7TJ, an d JZ3> lrom Appendix, we have for i ^ j 

(B ix (0)B jx (r)) = ~a 2 /12. (40) 

All B-correlation functions for time lags greater than r are zero. To understand this, consider (Bi X (Q)Bj X (2T)) . 
Averaging over 8^ T , 

(B ix (0)B jx (2T)) = ([A jx (r) - (m + l/2)a + sb, (2r)] [A£ fa + A&(-t) - rv ix ]), (41) 

where ma is the cell coordinate of particle j at t = 2r. The second term in (|41|) has no dependence on S2 T , while 
average of the first term gives zero. Again, this requires that the probability of any given configuration in a shifted 
cell is independent of ^r- 

These results can be summarized by the relation 

a 2 

(B la (nT)Bj/3(mT)) = — 5 al3 (l + S, tj ) [28 n>m - <5 n , m+ i - <5„, m -i] . (42) 

Fig. n presents simulation results for (B iy (Q)B iy (t)) in d = 2 for various collision angles a and a range of mean- free 
paths. Fig. contains corresponding results for (Bi y (0)Bj y (t)) . In both cases, the agreement with result (|42JI is 
excellent. Simulation results for (vi X (0)vi X (t)Bi y (0)Bi y (t)) as a function of time are presented in Fig. [3] for a similar 
range of parameters; the results are in agreement with the prediction of subsection IV Al that this autocorrelation 
function factorizes, and that the resulting ^-correlations are given by Eq. (|42|l . 

It follows that there are only two — a pure kinetic and a pure rotational (or collision) — contributions to the transport 
coefficients. Relation (|42|) is of central importance, because it contains all the geometrical features of the grid that 
contribute to the transport coefficients, and is independant of specific collision rules and particle properties. Since 
the kinetic contribution to the stress tensor is symmetric and has been calculated elsewhere, we concentrate here of 
the (collisional) contributions arising from B-correlations. 



C. Viscosities 



Explicit expressions for the collisional contributions to the viscous transport coefficients can be obtained by consid- 
ering various choices for k and a and (3 in Eq. (|25[) and using (|3J). Taking k in the y-direction and a = (3 — 1 yields 

1 oc ' 

v rot = v rct + v rot = __ £ £(^0)^(0)^)^^)). (43) 

n— i,j 

Eq. (|43|l is the expression used in Ref. to determine the collisional contribution to the shear viscosity. v\ ot and 
i/™* are the viscosities associated with the symmetric and the anti-symmetric contributions to the matrix of viscous 
transport coefficients. 

Other choices for k and a and (3 yield expressions for other linear combinations of the transport coefficients. In 
particular, the choice k = (1, 0, 0) and a = /3 = 1 yields 

(1 + (d- 2)/rf)< + 7 rot + e rot = ^^-^ E 5>»(0)iM0M*(*)5«(r>T)>. (44) 

n— i,j 

However, because of (|42|l . the right hand side of (|44|) is equal to (|43|l . so that 

(1 + (d - 2)/dK ot + 7 rot + e rot = v rot - (45) 

Finally, for k = (1, 1, 0)/V2 and a = 1, /3 = 2, one has 

((d-2)/dK ot -^ ot +7'' oi = 0, (46) 

since the resulting stress-stress correlation functions are zero. These results imply that e rot — 0, and that the 
longitudinal component of 0, which is the viscous contribution to the sound attenuation, is v rot . Finally, using these 
results in jnj), it follows that the collision contribution to the macroscopic stress tensor is 

a™* = K ot + vZ*)df,v a = v rot dpv a (47) 
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up to a tensor G with vanishing divergence, dpG a p = 0, which therefore will not appear in the linearized hydrodynamic 
equations. The collisional contribution to the effective shear viscosity is therefore i/ rot , and the viscous contribution 
to the sound attenuation is also u rot , instead of the standard result, 2{d — \)v jd + 7, for simple isotropic fluids. The 
corresponding hydrodynamic equation for the momentum density is therefore 

= -Vp + p(y km + !/ roi )Av + ^-^ fem v(V • v), (48) 

where we have used the fact that the kinetic contribution to the microscopic stress tensor, a km in (|16J) . is symmetric, 
and 7 fcm =00. Note that there is no collisional contribution to the last term in Eq. I|48|) . For d = 2, the viscous 
contribution to the sound attenuation coefficient in a simple liquid is v + 7. The results of the current calculation 
are consistent with this result, with v = v kvn + v rot and 7 = (as expected for a fluid with an ideal gas equation 
of state). In d = 3, while the shear viscosity is still given by v = v n + v rot , the viscous contribution to the sound 
attenuation coefficient is Av km /Z + i/ rot , instead of 4(z/ fcm + z^ rot )/3 + 7. The sound attenuation coefficient of the SRD 
model in three dimensions is therefore slightly smaller than in a simple liquid. Note, however, that the hydrodynamic 
equations are not affected. Although only certain linear combinations of the collisional contributions to the viscous 
transport coefficients are determined by Eqs. (f43|> — (f4fj|) . all coefficients in the linearized hydrodynamic equations are 
uniquely determined. 

Equilibrium measurements of the time-dependent density correlations jlGl ) in two dimensions yield results for the 
sound attenuation coefficient which are in good agreement with the theoretical predictions. A detailed comparison of 
these simulation results with theory will be presented elsewhere |l6j |. Currently, there are no similar measurements 
in three dimensions. 

The fact that the entropy of a fluid increases as a result of irreversible processes leads to certain positivity conditions 
on the transport coefficients [27| ■ SRD obeys an H-theorem p], Q , which implies that the entropy production is always 
non- negative. In Appendix E, it is shown, using a generalization of an argument from |27| . that the requirement of 
a positive entropy production leads to the conditions v km + v\ ot > 0, v\ ot > 0, and 7 > 0. Note that the result of 
Pooley and Yeomans [ijj for the collisional stress tensor a r °^ — v rot dpv a amounts to assuming v\ ot — v™ 1 = v rot /2 
and 7 = v rot /d. 



dv 
~dt 



+ (v • V)v 



VI. EXPLICIT EXPRESSIONS FOR THE COLLISIONAL CONTRIBUTIONS TO THE TRANSPORT 

COEFFICIENTS 



A. Viscosities 



Using the results of the previous sections, the collisional contribution to the viscosity can be written as 
2tN k B T 



i ^ "i - o.A T E {(^(0)^(0))(B il) (0)g J -„(0))+2(^(0)^(r))(B il) (0)B 3 - i< (r))} . (49) 



It is straightforward to evaluate the various contributions to the right hand side of 1)490. In particular, note that 
since velocity correlation functions only at equal time and for a time lag of one time step are required, molecular chaos 
can be assumed when evaluating these contributions, since it was shown in |Tlj that additional correlation effects only 
occur for larger time lags. Using (|42|l . the first term on the right hand side of Eq. 149fl reduces to 



,2 n 2 



Y, (Biy(0)B jv (0))(v ix (0)v jx (0)) =y^(i(0)) +jEE(«b(Ofe(0)) • (50) 
Momentum conservation, namely ^2 k Vkx(0) — 0, can be used to write 

EE(»«(%'(°)) = -E< u -(o)>, ( 51 ) 

i j^i i 

so that the right hand side of (|5(J|) reduces to 

|E^(°)) = i ¥BT ' (52) 
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where Nf = N — 1, because of momentum conservation. Similarly, using momentum conservation at t = r, the second 
term in on the right hand side of H49fl reduces to 

N 2 

J2 (B ty (0)B jy (r))(v lx (0)v JX (T)) = J>«(Q)t;i»(T)) • ( 53 ) 

1,3=1 i 



v rot _ » J j _ (^x(O)^(T)) . + Q( 



It follows that 

12r t fc B T 

For d = 2, if there are particles in the collision cell £f (t), the ensemble average of the term in brackets in H54fl is 

l-T-^<Wfa(0)vi«(T))U =(l-l/mi)[l-cos(a)], (55) 

where a is the collision angle. In three dimensions, the corresponding expression is 

1 2 

1 - 7^(%(0)%W)|mi = o( X ~ V"^)! 1 - cos(a)] , (56) 

for both models A and B. To obtain the final result, we now need to average over the number of particles, m.j, in the 
collision cell. If the average number of particles per cell is M, the probability that there are rrii particles in cell £f 
is given by the Poisson distribution P p (rrii,M) — e~ M M mi /mil. The corresponding (normalized) probability that a 
given particle, i, is in a cell containing a total number of particles is m,iP p {m,i, M) / M . Averaging now over the 
number of particles in a cell, we have, finally, 



a 2 (M-l + e- M 



Mt V M 



[1 - cos(a)] , (57) 



for all the collision models we considered (the standard model in d = 2 and both models A and B in 3d). Eq. I|57|) 
agrees with the result of Refs. ^| an d ^3 obtained using a completely different non-equilibrium approach in shear 
flow. Simulation results for the rotational contribution to the viscosity, Eq. I|57|) . are compared with the theoretical 
prediction in Fig. |l]for small M = 3, where our earlier approximation pj for M 3> 1 would not be accurate. The new 
expression correctly describes the limit M — > 0, where the collisional viscosity should vanish. 

B. Thermal diffusivity 

The collisional contribution to the thermal diffusivity can be calculated in a similar fashion. In particular, taking 
k = x in (T3j! and (O, and using we have D T = D^ n + D r T ot with 

oo ' N 

D k T m = - n \ bT2 E E (K 2 (0)/2 - k B T] [v](nr)/2 - k B T] ^(Q)^(nr)) (58) 

P n—Q — l 

i N 

Or' = g Nk y2 E {K 2 (0)" J 2 (0))(B K (0)B, I (0)) + 2(« l 2 (0) ? , 3 2 (r))(B, I (0)B 3I (r))} . (59) 

p B 1,3=1 

The kinetic contributions to the thermal diffusivity were calculated previously in both 2d and 3d ^Jj . Using the 
results presented earlier in this paper, it is straightforward to evaluate the collisional contribution to the thermal 
diffusivity. Just as momentum conservation was used to simplify the calculation of the collisional contribution to the 
viscosity, energy conservation, J2k=i v t = const., and ij^jl can be used to show that (|59^l reduces to 

Dt = 2^\ 1 -d(d+2)(k B Trj ■ (60) 

For d — 2, if there are mi particles in the collision cell £f (t), the ensemble average of the term in brackets in l|tj(J|) 



8(fc B T)2 



I 













(61) 
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In three dimensions, the corresponding expression is 

1 K 2 (oK(t)> 



15(fc B T) 2 



15to; 



— [1 - cos(a)] , 
to,- 



(62) 



for both models A and B. 

Using these results in Ij60|l and averaging over the number of particles in a cell, assuming again that the probability 
of having a given particle, i, is in a cell containing a total number of particles to.; is to.;P p (to.;, M)/M, where M is the 
average number of particles per cell, one finds 



Dl, 



1 



1 - e- M 1 



-Mi 



3(d+2)rM 



large M 3(d + 2)r M 



. M - -M 2 
small m 12(d + 2)r V 9 



- 1 






dx^ 


a; 




1 


1 


M " 


~ M 2 



[1 — cos(a)] 



[1 — cos(a)] 



[1 — cos(a)] 



(63) 
(64) 
(65) 



for all models considered. Note that in contrast to the viscosity, the rotational contribution to the thermal diffusivity 
is 0(1/M) for large M. Simulation results for the collisional contribution to thermal diffusivity are compared with 
(|63[) in Fig. |SJ This contribution to the thermal diffusivity, which is not negligible for small M (such as M — 3 in Fig. 
[5J, was not discussed in Refs. 01 an d E3- 

A comparison of the relative size of the rotational and kinetic contributions to Dt is given in Fig. H3 Since 
is independent of temperature while D!p n increases linearly with temperature, there is a temperature (or mean free 
path A = Ty/ksT) at which both contributions are equal. The ratio A/a of this specific mean free path to the cell size 
a is plotted as a function of the rotation angle a. Results obtained using Eq. H63|) . Eq. (89) of and Eq. (48) in 
[To| are presented for M = 5 in both two (dashed line) and three dimensions (solid line). In two dimensions, D™* can 
be larger than D!p n already at mean free paths as large as 0.25 a for large a. In d = 3, the rotational contribution is 
slightly less important. 



VII. CONCLUSION 

It has been shown that the random shift procedure introduced in Refs. 0,0 not only restores Galilean invariance, 
but also enables an exact evaluation of the collisional contribution to the transport coefficients. The current approach 
justifies in detail several assumptions used in the non-equilibrium calculations of Ref. jl2^ and [Tflj . and was used to 
determine the collisional contribution to the shear viscosity, the bulk viscosity, and the thermal diffusivity. 

A detailed analysis of the consequences of the fact that SRD collisions do not conserve angular momentum was 
also presented. It was shown that while the long-time, long-length-scale hydrodynamics of the model are not affected, 
it does lead to small changes in the viscous contribution to the sound attenuation coefficient. Although it has been 
pointed out previously 0] that the collisional contribution to the macroscopic viscous stress tensor is not symmetric, 
our interpretation of the consequences of this fact is different from that of Ref. 0|. In particular, the resulting slight 
modification of the coefficient of sound attenuation has no consequences for most practical applications, such as those 
in Refs. [lJll23; an d does not restrict the validity of the model. 
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IX. APPENDIX 

A. <(A£g 2 > 

If X s = Xi + S, with < Xi < a, 



r a/2 
%/2 

and 



(AQ X \ X ) 5 = [Q(-X s ) - G(X S - a)]dS (66) 

J -a/2 



Integrating over X, we have 



((Atlflxh = a f' 2 [Q{-X s ) + 0(X S - a)]dS. (67) 

J -a/2 



ra/2 r-a r-0 t'~& 

ts \2\ / dS I ' ,r " ' JX Jr ~ 2 

JO Ja-< 



((AZIY) = I dS I dx t + I d6 dx t = a z /A. (68) 

-a/2 JO 



B. i/j 



= i J a ' 2 ^ d8 dxi[Q(- Xi -6)- Q( Xi + S - a)] J . 



r a/2 ( r a 
a/2 

The integral over Xi is 



f dxi[Q(- Xi -5)- Q{ Xi + S - a)} = 50(6) - 6G(-6), 
Jo 



so that 

r a/2 



(69) 



(70) 



= -( d5[5*B(5) + 8 2 Q(-6)} = a 2 /12. (71) 

O J -a/2 

C. Various other correlations 

There are a number of other useful relations which are required to evaluate the B-correlations which can be easily 
evaluated using the same techniques. They include 



and 



and 



for i ^ j. Finally, 



{A8 x v ix ) = (A&{-T)v ix )=0 (72) 

(A&xA&s) = (A£i x A£j x ) = 0, (73) 

(*% x v jx ) = (Aa c (-T)v jx ) = 0, (74) 

<A&A&(-r)) = <A£ te A&(-r)) = (75) 

(A^(-r)[Ae fa (0) - A&(0)]) = -a 2 /12. (76) 
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To show this, assume < Xi(0) < a and na < Xi(0) + TVi X (0) < (n + l)a. Averaging first over S at time r and then 
the random shift at time t = 0, for fixed Xi(0) and 1^(0), one finds 

(A&(-T)[A£ fa (0)- A£? x (0)]| Xi(0)l x) = ((a/2-x l (0))(-a/2 + x l (Q)-Tv lx (0)). (77) 

The final ensemble average in 177|) reduces to 

I l-a 00 Mn+l)a-x\/T 

— I dx (a/2 — x) I dv x [— a/2 + x + tv x ]w(v x ) 

a JO n= -oo J(na-x)/T 







dx(a/2 — x) / dv x [~ a/2 + x + tv x ]w(v x ), (78) 



where we have dropped the index i and the time argument of x and v x for brevity. The integral over v x can be 
performed immediately, and the remaining integral over x gives the result, — a 2 /12. 

D. Proof of relation (vj X (T)vjy(T)v ix (0)B iy (0)} — 0. 

In order to evaluate this expression, averages over the random shift at time t — 0, <5o = {^Ox, $0y)i the shift at time 
t = t, S T = (5 TX ,S Ty ), and over the initial positions, i*j(0) = (xi(0), ?/i(0)), and velocities, Vj(0) of all particles are 
required. Averaging first over <5 , keeping all the other quantities fixed, yields Eq. i|3f])l. Next, note that £L(t) has 
an implicit dependence on the initial positions and velocities at t — and S T . We therefore write 

^(r)=^(r;<y r) {r fc } > {v fc ». (79) 

Because of translational symmetry, 

^(r; 6 r , {r fe }, {v fe }) = ^(r; 0, {r fc }, {v fc }), (80) 

with = Yk + 5 T . Keeping {r^} fixed, the average over 5 T in i|3U|) then become 
N 



iv poo poo r a /2 

"II / f({vk})dv k / dv k / Vj x (T)vjy(r)vi X (0) [£i V (T) - y 4 (0) + a/2 - Tv iy (0)] d6 T 

a j c —lJ— co J — oo J — a/2 

- II / /({ v fc» dv k / rff fe / u jx (r)« iy (T)w te (0) [£f y (r) - y 4 (0) + <5 r + a/2 - ru^O)] d5 T 

a k _ 1 J-oc J-oo J— a/2 

I!/ fd v k})dv k dT k v jx (T)v jy (r)v ix (0)[^ y (r)-yi(0) + a/2-TV iy (0)] . (81) 



where / is the N-particle Boltzmann distribution. 

The remaining average over the initial configuration can be split up into a sum of several terms. Each term 
corresponds to a situation in which the particle labelled i is restricted at time zero to be in a specific cell £i together 
with k\ other particles with given labels, while particle j is likewise residing only in a given cell £2 together with a 
set of k2 distinguishable particles. These restrictions are needed in order to have the post-collisional velocities Vj(r) 
and the cell label £| (r) to be independent on the initial positions of the particles in every term. 

We will show now that all these terms will vanish independently. Keeping the initial velocities fixed, an average 
over j)i(0) is performed under the condition mentioned above, i.e. where Cj y ( T ) = £i- These constraints mean that 
yi(0) must be in the interval [£1 — TVi y (0), £1 + a — rv iy (0)]. Since the post-collisional velocities of particle j are not 
affected by the position ?/i(0), each of the terms are proportional to 



£l+a-t>i y (0)T 



£i-v iy (0)i 



" ' ,- ■ a j/i(0) 



2 



£l+a— Tv iy (Q) 

= 0. 



/({vi}) dvf / & + - - w(0) - rv iy (0) dUO) cx a[6 + - - rv iy (0)] 

1 J 1 z 6-t^ k (o) 

All terms vanish independently, so that {v j X (r)v j V (r)vi X (0)5^(0)) is zero. Note that the same argument applies when 
£2 = £1, i.e. the particles being in the same cell, for both j ^ i and j = i. 
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E. Entropy production 



The general equation of heat transfer is [27 



'/ ( |^+v-Vs) ,<),<-„ -:- VlKVr.i. 



(83) 



where the stress tensor a a p is given in Eq. (J5J, s is the entropy per unit mass, p is the mass density, and k is the 
thermal conductivity. It is easily verified that a a p dfjV a can be written as 



v a p dpva = y ( dpv a + d a v (i - ^<S Q/ 3<9 7 v 7 j + y (V x v) 2 + 7(V • v) 2 



(84) 



where all three terms are independent, invariant under coordinate transformations, and can vanish independently. 
Using H83fl and (|84|l . it can be shown J27| that the entropy production in a given volume f2 is 



- / P sdV = 
at Jn j n 



+ * U Va + daVp _ 2 -S a ,d,v,)\ ||(V x v) 2 + 1(V • v) 2 



dV. 



(85) 



The requirement that the entropy production is non-negative implies that all transport coefficients in (|85|l are greater 
than or equal to zero, namely, 

v\ > 0, i> 2 > 0, 7 > 0, and k > 0. (86) 
For i>2 — this condition is identical to the well-known result discussed in [27|. For SRD, conditions l|86(l reduce to 

i/**** + v\ ot > 0, i/J"* > 0, and 7 rot > 0. (87) 
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FIG. 1: (Biy(0)B iy (t)} as a function of time. Results for A/a = 0.05, 0.10, 0.50,1.00, and for collision angles a = 60°, 90°, and 
120° are plotted, indicating that there is no dependence on the value of the mean free path. Time averages over 10 6 iterations 
were used to obtain the data. Parameters: L/a = 32 and M = 5. 
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FIG. 3: {vi x (0)vix(t)Biy(0)Biy(t)) as a function of time. Identical results are obtained for A = 0.05, 0.10, 0.50, and 1.00. Solid, 
dotted and dashed lines correspond to collision angles a = 60°, 90°, and 120°, respectively. Time averages over 10 iterations 
were used to obtain the data. Parameters: L/a = 32 and M = 5. 




FIG. 4: The normalized collisional contribution to the kinematic viscosity, u rot r/a 2 , as a function of collision angle a. The 
solid line is the theoretical prediction 15711 . The data were obtained by time averaging over 360,000 iterations. Parameters: 
L/a = 16, A/a = 0.1, M = 3 and r = 1. 




FIG. 5: The normalized collisional contribution to the thermal diffusivity, D^r/a 2 , as a function of collision angle a. The 
solid line is the theoretical prediction l|63[l . The data were obtained by time averaging over 360,000 iterations. Parameters: 
L/a = 16, A/a = 0.1, M = 3 and r = 1. 
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a 

FIG. 6: The specific ratio, A/ a, at which D!p n = D r T ot as a function of the rotational angle a for M = 5 in d — 2 (dashed line) 
and d — 3 (solid line). A = r^ksT is the mean free path, a is the cell size. For mean free paths below the curve the rotational 
contribution to the thermal diffusivity, _D y f , is larger than the kinetic part, D!p n . Theoretical expressions given by Eq. I63L 
Eq. (89) of Ref. 0, and Eq. (48) of Ref. |Hj were used. 



